Spin-orbit coupling assisted by flexural phonons in graphene 
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We analyze the couplings between spins and phonons in graphene. We present a complete analysis 
of the possible couplings between spins and flexural, out of plane, vibrations. From tight-binding 
models we obtain analytical and numerical estimates of their strength. We show that dynamical 
effects, induced by quantum and thermal fluctuations, significantly enhance the spin-orbit gap. 



I. INTRODUCTION 

Graphene combines in an unique way structu ral end 
electronic properties not found in other materials^]. Its 
special electronic structure made graphene the initial 
model for a two dimensional topological insulator P^l Nu- 
merical calculations suggested that the gap which is the 
hallmark of a topological insulator is in the micro Kelvin 
range, too small to be experimentally observed?*^ Re- 
cently, different approaches based on heavy adatoms de- 
position ha ve bee n proposed in order to engineer a topo- 
logical stateP^ELj Experimentally, these proposals are dif- 
ficult to be implemented, since the presence of adatoms 
implies a z — > — z (mirror) symmetry breaking which in- 
duces a Rashba-like SO coupling!^ These calculations 
are also the main guide for the interpretation of spin 
transport experiments. The calculations assumed lattice 
made up of ions of infinite mass, and ignored the lattice 
degrees of freedom. 

In the present work, we analyze the coupling of spins 
to lattice vibrations in graphene, with emphasis on flex- 
ural modes, which are unique to a two dimensional mem- 
brane, and have lower frequencies than conventional in- 
plane phonons. The main results of our paper are: i) 
A complete characterization of all possible couplings to 
flexural modes, including analytical expressions for each 
of them, and ii) An analysis of the effect of dynamic flex- 
ural modes on the Kane-Mele spin-orbit coupling, which 
defines the range of parameters where the topological in- 
sulator features of graphene can be observed. This anal- 
ysis is relevant in order to study the effect of ripples in 
spin transport! 1 ^ Our analysis is also applicable to carbon 
nanotubesp^ an issue that has become more interesting 
since carbon nanotubes have been recently propo sed a s 
a possible platform for hosting Majorana fermions! 15 * 16 ! 

We first present a symmetry analysis of the possible 
couplings between spins and lattice vibrations. This 
study can also be applied to static out of plane defor- 
mations. Next, in Section III, we make numerical es- 
timates of the couplings, based on simple tight binding 
models. Section IV discusses the change induced in the 
spin-orbit couplings due to the presence of lattice excita- 
tions, due to quantum and thermal fluctuation. Finally, 



Section V contains a summary of the more relevant re- 
sults. A number of mathematical details are included in 
the Appendices. 

II. SYMMETRY ANALYSIS AND COUPLINGS 

The low energy sector of the electronic spectrum lies 
on the two inequivalent corners of the Brillouin zone K±, 
known as valleys or Dirac points. In our convention 
K± = ± (471-/(3^), 0), where a is the carbon-carbon 
distance, see Fig. [T] The electronic Hamiltonian reads: 

n = -iv F f, ■ 8 + A/S 2 (g) s z + A fl (S x <8> s y - T, y ® s x ) 

(1) 

where X = (S x , T, y ) and S z are 4x4 matrices associated 
to the sublattice degree of freedom. This Hamiltonian 
operate in a space of 8-component Bloch functions ^ — 

Wu,K+,t> V>B,K + ,t, ^B,K_,t) -lpA,I£-,t, VU,K + 4, i>B,K + ,l, 

ipB,K.-,i,—'4>A,'K-,i) T - The first term corresponds to 
the massless Dirac Hamiltonian, the second term is the 
Kane-Mele mass which arises due to the SO interaction, 
and the last one is a Rashba-like coupling which is 
present in the case of a mirror symmetry breaking. Both 
A/ and A^ are weak, 1 — 15 /xeV in the former case 
according to previous estimates! 6 * 8 * 9 ! 



We analyze first the SO-coupling assisted electron in- 
teraction with flexural phonons from a symmetry group 
theory approach.^ Before we present the exhaustive 
analysis, we disclose the results in the following lines. 
The spin-phonon coupling Hamiltonian reads: 

T~L s -ph = H Al + T~Lb 2 + %c (2) 

where: 

H Al = 9i (Sa ® s y - Y* v (g) s x ) did l u Al + 
+92 [-A a ® s y (d 2 x u Al - dyU Al ) + 2A Z ® s x d x d y u Al } + 
+93 [(Sa ® s y + S y ® s x ) (d 2 x u Al - dyU Al ) + 
+2 (Sj. ® s x - Y< y <E) s y ) d x d y u Al ] 

(3) 
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TABLE I: Classification of the possible SO coupling terms 
according to how they transform under the symmetry oper- 
ations of Cq v . Note that all these operators are even under 
time reversal operation t — > — t. 



Ub 2 = 94^z ® S z (ub 2 )' 2 



(4) 



Hg' = .95[(A X ® s x )wi + (A^ ® s y )u 2 - 

-(A a ® Sy)u 3 - (Ay <E) s x )u 4 ] (5) 

The main ingredients in order to construct these 
couplings, electronic operators and symmetry-adapted 
phonon fields, are summarized in Tab. [i] and Tab. [n] 
respectively. Next, we construct these couplings step by 
step. 



Electronic operators 

The point group of the graphene crystal is C§ v , which 
contains 12 elements: the identity, five rotations and six 
reflections in planes perpendicular to the crystal plane. 
Instead of dealing with degenerate states at two inequiv- 
alent points one can enlarge the unit cell in order to 
contain six atoms, in such a way that K± are mapped 
onto the T point (see Fig. [l]b)). From the point of 
view of the lattice symmetries, this means that the two 
elementary translations (i ai , i a2 ) are factorized out of 
the translation group and added to the point group C§ v , 
which becomes C$ v — Cq v + t ai x Cq v + t &2 x Cq v . If 
we do not consider the spin degree of freedom, the 7r 
electronic states at K± transforms according to the 4- 
dimensional G' irreducible representation of Cq v . Then, 
in order to write down the low-energy electronic Hamil- 
tonian we have to consider the 16 Hermitian operators 
acting in a 4-dimcnsional space. We define two different 
sets of 4 x 4 hermitian matrices {Si}, {A^} associated 
to sublattice and valley degrees of freedom respectively, 
so the set {I, £j,Ai,£i • Aj} provides a representation 
of the algebra of generators of 11(4). In the basis in- 
troduced before XL and A.; matrices are odd under the 
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TABLE II: Classification of flexural (Z) phonon modes accord- 
ing to how they transform under the symmetry operations of 
Cg v . The numbers indicate the out-of-plane displacement of 
the atom within the 6-atoms unit cell. 



time reversal operation t — > —t, see Appendix A. If we 
consider the spin degree of freedom then the tt electronic 
states at K± transform according to the 8-dimensional 
representation G x -D1/2, where -D1/2 is the spinor rep- 
resentation associated to the spinorial part of the wave 
function. Thus, we introduce a set of Pauli matrices {sj} 
associated to the spin degree of freedom, also odd under 
time reversal operation. Importantly, we are now intro- 
ducing a pseudovector in the 3-dimensional space, mean- 
ing that the operators which contain s z are even under a 
reflection in the graphene plane (mirror symmetry oper- 
ation, z — > —z), whereas the operators which contain the 
in-plane components are odd. The result of this sym- 
metry based approach is summarized in Tab. [TJ where 
all the possible time reversal symmetric terms which in- 
volve spin operators are classified according to how they 
transforms under the symmetry operations of Cq v . These 
matrices correspond to all the allowed SO coupling terms 
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FIG. 1: a) Real space lattice and Brillouin zone of flat 
graphene with 2 atoms per unit cell, b) The same with 6 
atoms per unit cell. Note that K± are now equivalent to T. 



within the low-energy description. 



to construct the couplings. These are: 

|1> = 5 [- \AK + ) + \AK_) \BK+) + \BK_)} 

12} = \ [\AK + ) + \AK_) - |BK+) - |BK_>] 
|3) = l -[- \AK + ) + \AK-) + \BK+) \BK_}] 

|4) = \ [|AK + ) + \AK_) + \BK+) + |BK_>] (8) 

The polarization vectors of Eq. (|6| , together with he vec- 
tors of Eq. ([8]), form a symmetry adapted basis, in such a 
way that the displacement vector of a flexural mode can 
be written as \h) = UA t \A±) +u B2 \B 2 ) +u± |1) + u 2 |2) + 
it 3 1 3) + U4 1 4), where itj are the symmetry adapted (real) 
displacement fields. The results of this analysis are sum- 
marized in Tab. |TT] 



Phonon modes 

In the lattice with 6 atoms per unit cell, the am- 
plitude of a flexural phonon mode at T (= K±) is 
given by a 6-component vector whose entries are asso- 
ciated to the displacements of each sublattice atoms: 
\h) = (hAi,h B i,hA2, tiB2,h,A3, h B z)- This vector belongs 
to a 6-dimensional representation of C§ v which can be 
reduced as A\ + B 2 + G' . The polarization vectors as- 
sociated to the 1-dimensional irreducible representations 
correspond to the acoustic (ZA) and optical (ZO) modes 
at the original T point (with 2 atoms per unit cell) , whose 
polarization vectors in the 6 atoms basis read: 

i^i> = -4(i,i,i.M,i) 

V6 

\B 2 ) = -4(1,-1. 1,-1, 1,-1) (6) 
v 6 

The 4-dimensional irreducible representation corresponds 
to the 4 degenerate modes at K± points, whose polariza- 
tion vectors read: 
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Note that \A/B,K-) = (\A/B,K+))*. However, we 
must consider the real linear combinations of the vec- 
tors of Eq. (Fn) which transforms according to G' in order 



Spin-phonon couplings 

This analysis allows us to identify the SO assisted elec- 
tron coupling with flexural phonons at the center and 
the corners of the Brillouin zone. Since a uniform trans- 
lation of the crystal cannot affect the electron motion, 
it is clear that it can couple to acoustic phonons (at T, 
A\ phonons) only through spatial derivatives of the cor- 
responding displacement field. Moreover, if we consider 
the graphene as a continuum surface, only out-of-plane 
distortions which generate extrinsic curvature can couple 
to electron spin through the SO interaction. This implies 
that the leading term must depend on second derivatives 
didjUA 1 - Since did l UA 1 transforms according to A\, and 
(d^UAx — dyUA 1 , 2d x d y UA 1 ) forms a doublet which trans- 
forms according to E 2 , in principle three different cou- 
plings are allowed by the symmetries. We obtain Eq. 

m 

In the case of optical phonons at the center of the Bril- 
louin zone, note that there is no term which transforms 
according to the B 2 irreducible representation in the col- 
umn of z — > — z asymmetric operators, which means that 
the leading term must be quadratic on the phonon dis- 
placement fields. Since B 2 x B 2 = A\ we obtain a Kanc- 
Mele-like coupling term, Eq. |4j 

Finally, in the case of flexural phonons at the corners 
of the Brillouin zone, where both acoustic and optical 
branches are degenerate, the coupling reads as Eq. [5] 

III. TIGHT-BINDING MODEL 

The strength of these couplings can be estimated from 
a tight-binding model. We assume the convention of Fig 
[T] We choose the simplest tight-binding model with 4 
orbitals {s,p x ,Py,Pz} per carbon atom and only nearest 
neighbors hoppings. We neglect the effect of d orbitals 9 
because all these couplings are first order in the SO inter- 
action constant, as we are going to see. The Hamiltonian 
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TABLE III: Two-center matrix elements in the Slater-Koster 
approximation as function of the tight-binding parameters 
(see the text), p; represents a unitary vector in the direc- 
tion of maximum amplitude of the orbital pi, and 8 is the 
vector which connects neighboring sites A and B. 



can be written as: 
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(9) 



where the latin indices label the sites of the carbon atoms 
and A = s,p x ,p y ,p z labels the orbitals considered in the 
calculation. For the on-site energies we take t s = e s , 
t 



Pi 



e p . The two-center matrix elements can be com- 
puted within the Slat er-Koster approxim atio n as it is in- 
" We take (in eV)P 1121 e s = -7.3, 



III 



dicated in Tab 
e p = 0, V ssrT = -3.63 



ppll 



= -2.24. 



and Vpp a = 5.38, and 
We describe the spin-obit interaction in 



V — 4 2 

v spa L ±.^. 



terms of the Hamiltonian Tiso — A.L ■ s, where L and 
s correspond to the orbital angular momentum and spin 
operators respectively. We take A = 20 meVP" 



Our aim is to estimate the strength of the effective 
couplings within the low energy sector of the electronic 
spectrum. We can define the tt (p 2 , low energy sector) 
and a (s, p x , p v , high energy sector) orbital subspaces at 
least locally (see the discussion next and Fig. |2j. Then, 
the electronic Hamiltonian can be written in the block 
form: 



n 



^7T ^7T£7 



(10) 



We project out a orbitals by a Schrieffer-Wolf 
transformation:^ We take the Green function Q = 
(e — 'H) 1 , evaluate the block Q v associated to the low- 
energy sector, and use it in order to identify the low- 
energy effective Hamiltonian. If wc define Qf} a — 
(e — , then we can write: 



Git Q-na 
Qa-K Qa 



(11) 



5fl at : flat graphene 



S curv : bent graphene 




FIG. 2: Flat graphene and bent graphene. The arrows repre- 
sent the p x (black), p y (red), and p z (blue) orbitals. The local 
basis of atomic p orbitals in the bent graphene is uniquely de- 
fined by the isomorphism between them. 



We obtain Q v = 



(^ 0) ) 



so e 



On 1 = T-Ltt + H-KcrG^Ha-K- In the low energy sector (e 
0) the effective Hamiltonian reads: 







(12) 



Both the SO interaction and the out-of-plane distor- 
tions enter in the tt — a mixing blocks. In the ab- 
sence of out-of-plane distortions the only coupling al- 
lowed by the symmetries has the structure of a Kanc- 
Mele mass. Our estimation from the tight-binding model 
is A f / at = e s A 2 / (l8V s 2 pCT ) in agreement with Ref. HI We 
have a{'°* « 9 /ieV. 

In the case of B2 phonons we have to add to this anal- 
ysis the effect of a vertical displacement of one sublattice 
respect to the other. If the calculation is performed with 
the 6 atoms unit cell one can identify both the coupling 
with with flexural phonons at K± (C). In the case of 
A\ phonons, since the coupling depends on the second 
derivatives of the phonon field, the calculation is not so 
straightforward. 



Phonons at T 

As we mentioned before, the coupling with acoustic 
phonons can be inferred from the effect of extrinsic cur- 
vature of the graphene sample. The approach that we 
present here is quite similar to the calculation of the SO 
coupling in carbon nanotubes. The crucial fact is how 
to choose the basis of it and a orbitals. In the carbon 
nanotubes calculation the a orbitals are chosen in such a 
way that they follow the shape of the nanotube, and the 
7T orbital is chosen in the radial direction. Here we do es- 
sentially the same, but and we try to formalize it a little 
bit. In the low energy limit the graphene sample can be 
described within a continuum theory. We assume that 
the position of the carbon atoms lie on a smooth surface, 
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which is valid in the long-wavelength limit. Moreover, 
we assume that this surface is isometric to a plane. This 
means that a diffcomorphism / from a flat graphene sur- 
face (Sfi a t) to a curved graphene surface (S curv ) exists 
in such a way that the metric on the curved surface is 
pull-backed to the flat one. Thus, the isomorphism / de- 
fines an unique way to introduce a local basis for the p x 
and p y atomic orbitals. 

Consider graphene in a flat configuration. We intro- 
duce unitary vectors p x , p v in the direction of maximum 
amplitude of the orbitals p x , p y respectively, see Fig 
From a geometrical point of view, these vectors are e 
ements of the tangent bundle associated to <S/; a i- At 
the same time, the axis of maximum localization of the 
p z orbital verifies p z = p x x p y . Now consider a curved 
graphene which is related to the flat graphene by a iso- 
morphism /. The push- forward of / maps the tangent 
bundle of Sfi a t to the tangent bundle of S curv , which 
means that the vectors p x , p y at any position of the bent 
graphene are uniquely determined by the action of the 
push- forward of the isomorphism on the original p x , p y 
defined in the flat configuration. And of course, p z in the 
bent graphene surface is given by the vectorial product 
of the new p x , p y . More physically, what we are doing 
is the following. Consider a sheet of paper. Consider a 
point on it and write down a small arrow. Now deform 
the sheet of paper smoothly (folding is not differentiable) 
without breaking it. The new surface of your paper is 
then isometric to the original plane. Consider the previ- 
ous point in the new curved sheet of paper. The vector 
tangent to the arrow at that point, which is uniquely de- 
fined, is actually the push-forwarded original vector on 
the flat sheet of paper. This way of introducing the local 
basis in the bent graphene has two advantages: 1) we re- 
cover "smoothly" the original basis when we restore the 
curved graphene to the original flat configuration; 2) we 
keep the notion of parallelism by imposing / to be an 
isomorphism instead of just a diffeomorphism. Note that 
this apparent restriction does not affect the estimation of 
the SO assisted electron coupling with flexural phonons, 
since we are introducing extrinsic curvature and setting 
the Gaussian curvature to zero. 

We know how to choose the orbital basis and then cal- 
culate the two-center matrix elements, at least locally, 
using the Koster-Slater parametrization. For simplicity, 
we consider a curved graphene surface with a constant 
curvature along a given direction (a cylinder), so essen- 
tially the same problem as a carbon nanotube. Thus, we 
have two parameters, the radius of curvature R and the 
angle <f> between the direction of curvature and the x-axis 
(essentially the chiral angle in a nanotube), see Fig. [3] a). 
The new hoppings between 7r and a orbitals can be cal- 
culated following the prescription of Tab. |I1I| These are 
function of the angle 9 defined in Fig. [3]b). Assuming 
that R 3> a, we have to the leading term in a/R: 



a) 





e 



FIG. 3: a) Definition of the angle <f). b) Sketch for the calcu- 
lation of the new hoppings between p z and Pi orbitals. 



where 8 x ^ y are the components of the vector 6 which con- 
nects nearest neighbors. After a straightforward calcula- 
tion, the block Hamiltonian that mixes it and a states at 
K± can be written as the matrix: 
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where V\ = V, 



ppa 



V, 



ppTT 1 



(14) 

V 2 = {V PP a + 3V ppn ) /2, and 
t = ±1 labels the valley K±. 

This expression is exact to the leading order in a/R as- 
suming a constant R and (f> along the graphene surface. 
Now we perform a local approximation, which is valid 
at long wavelengths: we assume that R and <f> depends 
slightly on the position. Hence, they can be related with 
the second derivatives of the height profile, since the sec- 
ond fundamental form (J 7 ) of a surface in the Monge's 
parametrization (r = (x,y,h(x,y))) reads: 



T = 



1 

VI + dih&h 



d x d x h d x d y h 
d x d y h d y d y h 



(15) 



If we neglect quadratic terms on h, then J- is just the 
tensor of second derivatives. At the same time, note that 
within the continuum description of graphene as a mem- 
brane, the height profile h should be identified with the 
flexural acoustic phonon field at long wavelengths ua 1 ■ 
Then, the local R and <fi can be related with the second 
derivatives of Ua 1 as: 



d x d x u Al 



-iT 1 , 



dyd y u Al w -R 



sin 



d x d y UA 1 ~ —R 1 sin <j> cos <f> 



(16) 



\5 x cos{<t>) 



■8y$m(4>)\ 



R 



(13) 



By projecting out the a electronic states as Eq. [12] 
indicates, we arrive at the electron-phonon coupling of 
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Eq. [3] with: 



ae s A (Vppv + Vppn) % 

91 = — « 3me v.A 

g 2 = aVpp * A « 4 meV • A 

2 (Vppa V'ppTr) 



(17) 



Importantly, within the present Koster-Slater approxi- 
mation we have to extend the tight-binding calculation 
to second-nearest neighbors in order to obtain a non-zero 
33 coupling. In that case we obtain: 



aA (3Vpj 



53 = 



where Vppi- and Vpph are new second-nearest neighbors 
hopping parameters. 

In order to estimate the coupling with optical phonons 
we have to consider the effect of a vertical displacement 
of one sublattice with respect to the other, similarly 
to the case of siliceneP^ Processes which involve only 
one phonon give a vanishing contribution as expected 
from symmetry considerations. We can repeat the same 
scheme as before by considering virtual processes medi- 
ated by two phonons (oc Hl l ^ x %~ 1 %l l ^ x , note that H a 
also contains the SO interaction). We identify the Kane- 
Mele like coupling with flexural optical phonons at T, 
whose strength reads: 



V ppa ) (y, 



(2) 



ppa 



v. 



(2) 



PP7T 



VppTT ) 



(18) 



,(2) 



54 



2e?A(I/ pi 



V m 



9aW s % a 



20 meV • A' 



(19) 



Phonons at K± 

We estimate now the coupling with phonons at the cor- 
ners of the Brillouin zone. We simplify the tight-binding 
model in order to treat the problem analytically. We 
are going to consider the model described in Ref. |2H 
and adapted in Ref. [55] in order to describe the acoustic 
phonon modes in graphite. We consider two parameters: 
V on , which is the on-site energy of the a orbitals, and 
Vhopi which is the hopping between a orbitals at nearest 
neighbors when the orbitals are maximally localized in 
the direction which links the two atoms, otherwise the 
hopping is taken to zero. These parameters can be esti- 
mated from the Slater-Koster parameters as: 



Vhop — 



v., 
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2\f2V S p„ 
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2V r 
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(20) 



This model was employed in Ref. [5] in order to esti- 
mate the SO coupling in graphene and carbon nanotubes. 
Within this model for flat graphene, the Kane-Mele cou- 
pling reads Aj lat — 3V OIl A 2 /(4V^ op ), in agreement with 



Ref. IH1 Note also that this estimation is numerically very 
close to the one of Ref. [U 

We perform the calculation in the unit cell with 
6 atoms. We can estimate the terms that mix 7r 
and a states by considering the vertical displace- 
ment of one lattice respect to the other, as we 
mentioned before. By projecting out the a orbitals 
we arrive to a 6x6 effective Hamiltonian for the 7r 
electronic states at the new T, which can be seen 
as matrix expressed in the monoelectronic basis 
(|IM l7 r), \TA 2 TT),\TA 3 n), \TB l7 r), \TB 2 n) , \TB 3 *)). 
In order to identify the effective Hamiltonian in 
the low energy sector we have to express this 
matrix in the monoelectronic basis associated to 
the lattice with 2 atoms per unit cell, let's say 
(\TAn), \TBtt), \K + Att), \K + Btt), |K_Att), |K_Btt)). 
Both basis are related by the unitary transformation: 
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(21) 



By doing so, we identify the strength of the coupling with 
phonons at the corner of the Brillouin zone: 
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2V2aVl p 



4 meV • A 1 
(22) 



IV. ENHANCEMENT OF KANE-MELE MASS 

One of the most interesting consequences of this anal- 
ysis is the effect of the coupling of Eq. [4] in the electronic 
spectrum. The contribution of flexural phonons to the 
Kane-Mele coupling can be written as: 



A P h = gi({uB 2 ) 2 ) 



(23) 



where the brackets express the thermal average over the 
entire Brillouin zone. 

The flexural optical (ZO) mode strictly at T transforms 
according to B 2 . A rough estimate consists on neglecting 
the contribution from the acoustic branch and the disper- 
sion of the optical mode. Since hu^° « 110 meV, tem- 
perature plays no role, see the red curve in Fig. [IJ How- 
ever, the zero-point motion contribution Ao = 2M ^ 4 zo ~ 

0.03 meV is non-negligible (here M is the mass of the 
carbon atom). 

This is a very crude approximation, since the iden- 
tification of UB2 with the ZO mode is strictly true at 
the r point. The entire Brillouin zone contributes to 
the average, so away from T both flexural acoustic (ZA) 
and optical branches enter. We can use the symmetry- 
adapted basis \A\) and \B 2 ) phonons as a basis in order 
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FIG. 4: Effective Kane-Mele mass induced by the coupling 
with flexural phonons. In red (lower curve) the estimation 
neglecting the acoustic branch and the dispersion of the op- 
tical one. In blue (upper curve) the calculation within the 
model described in Appendix B. 



to describe the polarizations of the v = ZA,ZO phonons: 



\v)=rf Al (q)|Ax)+^ 2 (q)|B 2 ) 



(24) 



Note that, because of time reversal symmetry, r\ (q) = 
[rj (— q)]*. Thus, we have: 



<(«*) a >«4 S IX M 



where: 



A? 



v 2— 1 



(25) 



(26) 



is the Fourier transform of the phonon displacement field 
in branch i>, and the brackets denote thermal average. 

We need a model in order to describe the deviations 
of the polarizations vectors and the frequencies of both 
branches in the entire Brillouin zone. In Appendix B we 
describe the simpler nearest-neighbor forces model that 
one can consider in order to describe the dynamics of 
flexural phonons. We compute the induced Kane-Mele 



mass from Eqs. (23 1 and (25) within this model. The 
results are shown in Fig. |4| (blue curve). Remarkably, 
the Kane-Mele gap induced by phonons 2A p h is of the 
order of 0.1 meV. 



V. CONCLUSIONS 

We have studied the possible SO-mediated electron in- 
teraction with flexural phonons allowed by the symme- 
tries of the lattice and estimated the strength of these 
couplings from a tight-binding model. This analysis can 
be used in order to study the SO coupling in carbon nan- 
otubes. 

We find that the quadratic coupling with the phonons 
at r, particularly with the ZO branch, renormalizes the 



Kane-Mele mass in a remarkable way. From our the- 
ory, we predict an enhancement of two orders of mag- 
nitude in comparison to previous estimations for flat 
graphcnc, putting this gap close to the present exper- 
imental limitiPfi]. Note that the frequency of flexural 
modes depends on the coupling to the substrate, if any, 
and can be tuned by applied strains^. On general 
grounds, it can be expected that a compressive strain 
will lower the frequency of these modes, enhancing the 
spin-orbit coupling. 

Our theory is also relevant for spin transport experi- 
ments, particularly in suspended samplesp^l where charge 
transport is ultimately limited by flexural phononsP^l 
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Appendix A: Choice of the spinor basis 

With this appendix we try to avoid possible confusions 
due to notation and choice of spinor basis in the main 
text, in particular when compared to Ref. |H 

When discussing the symmetries of the graphene lat- 
tice from the perspective of group theory it is common 
to introduce electronic operators without specifying their 
explicit expression in a certain basis, only the algebraic 
rules which determine how they transform under the 
symmetry operations of the enlarged point group Cq v . 
For instance, the two 4x4 matrices which transform ac- 
cording to the vector irreducible representation E\ are 
denoted by S x , Y. y . The matrix E z is defined in order 
to complete the Pauli-matrix algebra, and it transforms 
according to Ai- On the other hand, the matrices which 
transform according to the vector irreducible representa- 
tion E' 2 are denoted by A x , A y , and as before A z is de- 
fined in order to complete the Pauli-matrix algebra and 
transforms according to B±. 



A more physical construction would be the following. 
As we mentioned in the text, in order to write down the 
low-energy electronic Hamiltonian (neglecting the spin) 
we have to consider the 16 Hermitian operators acting 



Irrep 


z 


— > — z symmetric 


2 — > — z asymmetric 


Ai 


O z ®T z ® S z 


<J X ® T Z ® Sy — Oy ® S X 


A 2 




X ® T Z ® S X + Oy ® Sy 


B 2 


r z <g> s z 




E 1 




/ Oy ® S Z \ 

\ O x ® T z <g> S z ) 






f -0 Z ®T Z ® Sy N 

y O z ® T z (g) , 




E 2 




( X ®T Z ®Sy+Oy®S X \ 
\0 X ®T Z ® S X — Oy® Sy J 
( -T Z ® Sy \ 

V Tz ® 8 X J 


E{ 


( 


' —Oy ® Ty ® S Z 
v Oy ® T x ® S Z ^ 


) 




G 






^ — (Ty ® Ty ® Sa; ^ 
— CTy ® Ty ® Sy 
<Ty (g> Ta; ® Sy 

\ o-y ® r x ® Sj y 





and then the possible SO terms are constructed from 
products of a spin matrix s, with £j or Aj. 



Appendix B: Model for fiexural phonons 

We consider the simpler nearest-neighbor forces model 
where the elastic energy of the lattice can be written as: 



E = 



E 



I J2 Hb A + [ hBi - \ J2 Ha j 

(ij) J \ («) / 



(Bl) 



TABLE IV: Classification of the possible SO coupling terms 
according to how they transform under the symmetry opera- 
tions of the lattice and reflection in the graphene plane. 



Here a has units of energy and can be related 
with the bending rigidity of graphene in a continuum 
description!^ This model leads to the dynamical matrix: 



in a 4-dimensional space. We could define two differ- 
ent sets of 2 x 2 Pauli matrices {<7j}, {t;} associated to 
the physical sublattice and valley degrees of freedom re- 
spectively, so the set {Z, <7j, Tj, erj <g) r.,} provides a rep- 
resentation of the algebra of generators of E/(4). These 
operators act in a space of 4-component Bloch functions 
* = (ipA,JS.+ ,ipB,K+,ipA,K-,t: 4>b,k.^) T - Then, we intro- 
duce the Pauli matrices associated to spin, so we have to 
double the space of Bloch functions in order to contain 
the two spin projections. The possible SO couplings in 
this basis are summarized in Tab. IIVI In the first row we 
can recognize the Kane-Mele and Rashba-like couplings 
as discussed in Ref. 0J Note that in this basis the time 
reversal operation is implemented by the anti-unitary op- 
erator T — ir x ® SyJC, where K, denotes complex conju- 
gation. 

In the main text we employ a different basis, 
where the order of the projection of the Bloch func- 
tions at each sublattice in different valleys is in- 
verted, and also a minus sign is introduced: = 
(Vm,k+,'0b,k + ,V'b,k_,-V'a,k_) T - This basis is very 
convenient essentially because the notation is simplified. 
In this basis the operators £j, Aj are related with the 
matrices acting in subalattice and valley indices as: 



A., 



Oj 



3 J 

n 



(Al) 



The notation is simplified essentially because in the basis 
the time reversal operation is implemented by T = iYi y ® 
Ay Cg) SylC, in such a way that the three sets of matrices 
are odd under the action of T: 





► ^y^i — 


-Si 




-A Ay A* Ay = 


-A, 








Si 


► s y s iSy = 


-Si 



(A2) 




+ ^ -2/(q) 
-2/(q)* 3+^ 



(B2) 



where / (q) = ^ Q e lq ' 5a , and the sum is extended to 
nearest-neighbors. The frequencies read: 



U)± 



I 2a 
3Ma 2 



l/(q) 



±2|/(q) 



(B3) 



with polarization vectors |±) = ( |/(q)| ' ' ^ 

q = we have ui§° = uo+ = yQf^~, ^ A = w_ = 0. We 
have for q ~ F: 



U), 



8a 
Ma 2 

ZA 



8M 



8Ma~ 



(B4) 



Note that the dispersion relation of ZA phonons is 
quadratic, as expected. Moreover, we deduce the rela- 
tion a = 6y/3n. At K±, w|| 
same time, we obtain: 



ZA 

± 



,zo 



f,°/2. At the 



K°(q)| = 



1 + 



g/(q) 
|/(q)l 



(B5) 



where 3?/ denotes the real part of /. We compute the 



induced Kane-Mele mass from Eqs. (23 1 and (25) within 



this model. We set the value of a from uj^P . We obtain 
a = 8.5 eV (re ps 0.8 eV). The results are shown in Fig. 
[1] (blue curve). 
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